AD-764  711 

PREDICTING  THE  DAY-TO-DAY  VARIABILITY 
OF  THE  MID-LATITUDE  IONOSPHERE  FOR 
APPLICATION  TO  HF  PROPAGATION  PREDICTIONS 

Charles  M.  Rush,  et  al 

Air  Force  Cambridge  Research  Laboratories 
L.  G,  Hanscom  Field,  Massachusetts 


23.  May  1973 


J5  m  Royal  Rosd,  Springfield  Va.  22151 


AFCRL.TR-7M>335 
23  MAY  1973 

AIR  FORCE  SURVEYS  IN  GEOPHYSICS.  NO.  268 


IONOSPHERE  PHYSICS  LABORATORY  PROJECT  8666 


AIR  FORCE  CAMBRIDGE  RESEARCH  LABORATORIES 


L.  G.  HANSCOM  FIELD,  BEDFORD,  MASSACHUSETTS 


Predicting  the  Day-to-Day  Variability  of  the 
Mid-Latitude  Ionosphere  for  Application 
F  Propagation  Predict 


CHARLES  M.  RUSH 
JOSEPH  GIBBS* 


^ARCQN  Coiporo*)an,  WoWdold,  Massochu**tts 


Approved  for  public  r»lou*e;  distributism  url(n»lt»d. 


AIR  FORCE  SYSTEMS  COMMAND 
United  States  Air  Force 


A  :- 


y  ■  •  .  ■  •  ■  ;  •'  '  .  '■‘  •s'.  ■  ■  v  ir  •  . .  •  -  -  y'  ■ 

.*•  Ji- v  .>  '' U  >.;%  .  ■  ;v*.  ‘ •-* ' .  ■.  \  \  ■ 

.  »  .  ■*  •  *  *■-  '  *i  -  -  .**  *•  4*  -•  ,  »»*-*.  1-.*'  .  ■■  .• 


Unclassified 

Security  CUiil fleet  loo 


DOCUMENT  CONTROL  DATA  •  RAD 

(Security  classification  of  title,  body  of  abstract  and  indexing  annotation  must  be  tittered  uAe*  tke  overall  report  is  classified} 


n  omomATiMQ  ACTIVITY  (Corporate  eviAorJ 

Air  Force  Cambridge  Research  Laboratories  (LII) 

L.  G.  Hanscom  Field 

Bedford,  Massachusetts  01730 


la.  eeeottr  mcuaity  cu*mric*TK>w 

Unclassified 

nr"asouu  —————— 


I.  UtM*T  TITUI 

PREDICTING  THE  DAY-TO-DAY  VARIABILITY  OF  THE  MID-LATITUDE 
IONOSPHERE  FOR  APPLICATION  TO  HF  PROPAGATION  PREDICTIONS 


4.  ouciwtiv*  mote*  (Typt  of  npert  md  i»ciu.\tin iaU>) 

Scientific.  Interim. 

s.  ^UTHORiS)  { A-'.c  v**e,  i/tiUU,  fan  nmttj 

Charles  M.  Rush 

Joseph  Gibbs 

t.  AIROAY  DATt 

23  Mftv I.fava  . 

?a  total  ho.  or  mqio 

HO.  0^  RIM 

'  9  _ 

<&  cortr&t  G*  GHAUT  «0, 

6.  RROdfCT.  TMK,  WORK  VMlT  Mb*  86660103 

A  OOO  JUMIUMHT  681000 

®<*  QHlfilHATO***  *I»94T  MUNttlflW 

AFCRL  ~TR  -73-0335 

**>  0TH«Rii|f O^THCTU (A>q awfcwt 

AFSGNo.  268 

«.  W*TW»V)Ytq*  *TAYfM«M?  " . 

Approval  for  public  release;  distribution  unlimited. 

M.  «U#»VI*«MTANV  ROTW 

TECH,  OTHER 

Til  ftJTWUV - 

Air  Force  Cambridge  Research 
Laboratories  (LII) 

L.  G.  Hanscom  Field 

Bedford,  Massachusetts  01730 

Ifc./lMHjVeACT 

'  Hourly  values  of  ground-based  ionosonde  data  Lave  been  subjected  to  a 
detailed  statistical  analysis  in  order  to  determine  tme  day-to-day  variability 
displayed  by  various  ionospheric  regions.  It  is  found  that  the  standard  de¬ 
viation  about  the  observed  F2  region  monthly  median  critical  frequency  is  two 
to  three  times  that  observed  in  the  E  and  FI  regions  at  mid -latitudes.  The 
day-to-day  variability  of  foE  and  foFl  is  such  that  monthly  median  values  of 
those  parameters  can  be  used  to  represent  the  daily  variation  in  the  E  and  FI 
region.  Title  would  imply  that  HF  prediction  programs  which  rely  on  median 
data  for  specification  of  the  ionosphere  yield  results  that  can  be  used  in  day- 
to-day  operations  for  those  propagation  modes  that  are  controlled  by  the  E  andj 
FI  region.  As  the  F3  region  ia  highly  variable,  a  method  of  using  timely 
observations  of  foF2  to  predict  the  expected  values  of  the  F2  region  critical 
frequency  has  been  developed.  The  method  has  the  advantage  that  errors  in 
predicting  the  monthly  median  foF2  can  be  averted  and  the  predicted  values  ar^ 
generally  closer  to  the  actual  observations  than  the  observed  monthly  median. 
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Abstract 


Hourly  values  of  ground-based  ionosonde  data  have  been  subjected  to  a 
detailed  statistical  analysis  in  order  to  determine  the  day-to-day  variability  dis¬ 
played  by  various  ionospheric  regions.  It  is  found  that  the  standard  deviation 
about  the  observed  F2  region  monthly  median  critical  frequency  is  two  to  three 
times  that  observed  in  the  E  and  FI  regions  at  mid -latitudes.  The  day-to-day 
variability  of  foE  and  foFl  ia  suoh  that  monthly  median  values  of  those  parameters 
can  be  used  to  represent  the  daily  variation  in  the  E  and  FI  region.  This  would 
imply  that  HF  prediction  programs  which  rely  on  median  data  for  specification  of 
the  ionosphere  yield  results  that  can  be  used  in  day-to-day  operations  for  those 
propagation  modes  that  are  controlled  by  the  E  and  FI  region.  As  ■  he  F2  region 
ia  highly  variable,  u  method  of  using  timely  observations  of  foF2  tc  predict  the 
expected  values  of  the  F2  region  critical  frequency  has  been  developed.  The 
method  has  tho  advantage  that  errors  in  predicting  the  monthly  median  foF2  can  be 
averted  and  the  predicted  values  are  generally  closer  to  the  actual  observations 
than  the  Observed  monthly  median. 
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Predicting  the  Day- to- Day  Variability  of  the 
Mid-Laiilyde  Ionosphere  for  Application  to 
HF  Propagation  Predictions 


1.  INTRODUCTION 

Routine  ground-based  ionospheric  soundings  have  over  th«  years  provided  a 
great  deal  of  information  concerning  the  structure  of  the  ionosphere  and  its  vari¬ 
ability  in  both  time  and  space.  These  observations  have  been  w  jd  to  help  deter¬ 
mine  the  various  physical  mechanisms  tiist  control  the  ionospheric  behavior  and 
provide  the  basis  for  ascertaining  how  changes  in  the  ionosphere  affect  propagation 
of  electromagnetic  energy  over  a  wide  band  of  the  frequency  spectrum.  Perhaps 
the  moat  significant  and  well-known  aspect  of  the  ionosphere  ia  its  apparent  control 
over  high  frequency  (HF)  radio  propagation  conditions.  As  too  ionosphere  is  the 
medium  in  which  HF  signals  propagate*  changes  in  the  ionospheric  structure  load 
naturally  to  changes  in  HF  propagation  choracteristios, " 

Because  of  the  bnportonffie  of  HF  communication  in  «  variety  of  civilian  curl 
military  ■uses,  &  great  deal  of  effort  has  bean  spent  on  predicting  the  structure  of; 
the  ionosphere  tor  the  express  purpose  of  improving  HF  propagation  predict  ions. 
For  the  most  part,  HF  prediction  techniques  rely  upon  median  bases  for  the  speci¬ 
fication  ar-<!  prediction  «s?  the  ionosphere  and  are  thus  only  strictly  valid  for  median 
conditions,  There  is,  however,  an  increasing  amount  of  off  davoied  to  applying 
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median  predictions  in  normal  day-to-day  applications.  This  would  lead  to  meaning¬ 
ful  results  only  if  the  median  behavior  of  th?  ionosphere  provided  a  good  represen¬ 
tation  of  the  daily  variability. 

In  this  report,  the  variability  of  the  ionosphere  about  its  median  value  is  in¬ 
vestigated  using  data  that  are  representative  of  various  ionospheric  regions. 
Although  there  are  other  factors  that  affect,  and  in  certain  circumstances  limit, 

HF  propagation  (absorption,  terrestrial  and  man-made  noise,  for  example)  only  the 
variability  displayed  by  the  ionospheric  electron  density  is  studied.  This  has  been 
accomplished  by  making  use  of  the  hourly  critical  frequency  tabulations  of  the  nor¬ 
mal  E  (foE),  FI  (foFl)  and  F2  (£oF2)  regions  observed  at  a  number  of  mid-latitude 
locations.  An  attempt  has  also  been  made  to  relate  the  observed  daily  variability 
of  the  ionosphere  to  the  practicality  of  U3ing  median  predictions  in  day-to-day 
operations. 

It  will  be  seen  that  the  monthly  median  values  of  foE  and  foFl  provide  a  rela¬ 
tively  good  representation  of  the  daily  values  of  those  particular  parameters  at 
mid-latitudes.  In  the  F2  region,  on  the  other  hand,  the  daily  variability  of  foF2 
is,  at  certain  times,  quite  large  and  the  monthly  median  value  of  foF2  appears  to 
be  of  little  practical  use  in  daily  operations.  We  have,  therefore,  investigated 
methods  by  which  this  daily  variability  can  be  predicted  assuming  that  timely  ob¬ 
servations  of  foF2  are  available.  It  will  be  shown  how  much  improvement  hi  the 
prediction  of  foF2  can  be  obtained  using  immediate  past  values  of  that  parameter 
to  update  a  predicted  value.  Whether  the  amount  of  improvement  that  is  obtained 
justifies  the  cost  necessary  to  implement  such  a  technique  must  be  decided  with 
reference  to  the  specific  application. 

&  VARIABILITY  OF  THE  IONOSPHERE  ABOUT  THE  OBSERVED  MEDIAN 

In  order  to  obtain  ''  reliable  estimate  of  the  dally  variability  of  the  ionosphere 
about  the  observed  monthly  median,  hourly  values  of  the  critical  frequencies  of 
the  E,  FI  and  F2  regions  have  been  obtained  for  a  number  of  mid-latitude  iono¬ 
spheric  sounding  stations  that  were  operating  during  the  years  of  1050  and  1964. 
These  tone  periods  were  chosen  to  be  representative  of  the  extremes  in  solar 
activity  (1858  -  color  maximum;  1064  -  solar  minimum).  Since  these  were  years 
when  international  agreement  concerning  all  types  of  geophysical  observations 
were  in  effect  UGY  and  IQSY),  ihe  number  of  stations  supplying  observations  tend 
to  bo  more  numerous  than  in  other  yearn,  Data  from  over  150  ionosonde  stations; 
that  wore  operating  the  entire  year  of  1950  were  used  to  det  rmine  the  variability 
of  foF3  during  solar  maximum  conditions,  while  name  40  mid -latitude  station*  were 
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investigated  for  1964  to  determine  the  variability  of  foF2  daring  solar  minimum 
conditions.  Hourly  data  for  foE  and  foFl  were  used  for  sill  stations  north  of  45° 
geomagnetic  latitude  for  the  months  of  March,  June,  September  and  December  in 
both  1958  and  1964,  Only- the  results  that  applyto  the  mid -latitude  ionosphere  will 
be  discussed  here. 

For  each  of  the  ionospheric  parameters  and  stations  available,  the  observed 
monthly  median  for  each  hour  was  determined  and  the  individual  hourly  observa¬ 
tions  were  then  subtracted  from  the  monthly  medians.  The  individual  deviations 
were  also  norm  alio  ad  to  their  respective  hourly  median  values  to  enable  study  of 
both  the  absolute  as  well  as  relative  variability  in  the  various  regions,  of  the  iono¬ 
sphere,  These  deviations  (both  relative  and  absolute)  were  then  subjected  to  a  de¬ 
tailed  statistical  analvsia  and  the  first  four  statistical  moments  were  computed  for 
each  data  set  (Kenney  and  Keeping.  1954).  A  standard  computer  program  was  used 
♦o  determine  the  mean,  standard  deviation,  skewness,  and  kurtosia  for  each  set  of 
deviations.  It  is  the  standard  deviation  that  provides  the  most  straightforward  in¬ 
formation  concerning  the  variability  of  the  hourly  observations  about  the  observed 
monthly  median.  In  order  to  increase  the  significance  of  the  results,  the  devia¬ 
tions  for  each  station -month  wave  grouped  into  three  hour -time  blocks  centered  on 
0000,  0300,  ...  3100  h  before  performing  the  statistical  analysis.  The  groupings 
were  done  both  in  local  mean  time  (LMT)  and  Greenwich  moan  time  (GMT), 

The  results  of  the  analysis  indicated  that  the  hourly  deviations  of  ioE  and  foFl 
could  reasonably  be  considered  as  a  normal  ensemble,  The  average  deviation 
about  the  median  rarely  exceeded  two  percent  implying  that  the  monthly  median  and 
the  monthly  mean  were  nearly  equal.  The  skewness,  which  provides  an  indication 
of  the  asymmetry  of  the  data  acts  was  also  rather  close  to  aero  as  it  must  be  for  a 
normal  distribution.  In  the  F2  region,  there  were  occasions  waer.  the  average  de¬ 
viation  about  the  median  was  on  the  order  of  ton  percent  and  the  deviation  data 
appeared  to  be  quite  skewed,  For  the  most  part,  this  was  found  to  be  associated 
with  periods  of  severe  magnetic  disturbances.  However,  to  a  first  approximation, 
it  is  felt  that  the  F£  region  deviations  could  also  be  considered  as  normally  dis  ¬ 
tributed  for  the  majority  of  eases  studied. 

The  statistical  analysis  indicates  that  the  variability  of  the  ionosphere  le  not 
only  different  for  the  different  ionospheric  layers  but  that  the  variability  in  any 
given  layer  displays  local  time,  seasonal  and  aolar  cycle  dependencies  of  varying 
degrees.  Figure  la  shows  the  standard  deviation  in  MH*  about  the  monthly  median 
of  ioE  tliat  U»  typically  observed  at  mid -latitude  location?.  Shown  In  the  figure  are 
results  obtained  from  data  taken  at  Ottawa,  Canada,  and  Slough,  England,  during 
March,  Jane  and  September  iS)58  and  1964,  The  geographic  and  geomagnetic 
coordinates  of  these  and  ether  stations  whose  results  are  shown  in  the  import  are 
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given  in  Table  1.  Because  the  E  region  is  under  strict  solar  control,  routine  ob¬ 
servations  of  foE  are  only  available  during  daylight  hours.  For  the  local  mean 
times  shown  (0900  to  1500  h),  the  standard  deviation  of  foE  about  the  median  is 
never  greater  than  0. 05  MHz,  with  the  difference  between  solar  maximum  and  solar 
minimum  results  being  small.  This  difference  is  greatest  during  March  and 
September  and  practically  nil  during  the  summer  solstice.  Although  not  shown  in 
the  figure,  the  standard  deviation  at  times  surrounding  0600  h  and  1800  h  tends  to 
be  slightly  larger  than  for  midday.  However,  the  scarcity  of  reliable  data  at  these 
times  renders  it  difficult  to  determine  how  statistically  significant  the  r  esults  are 
during  periods  surrounding  sunrise  and  sunset. 
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Figure  la.  Absolute  Standard  Deviation  About  the  Observed  Monthly 
Median  foE  at  Ottawa  and  Slough  During  March,  June  and  September, 
.1958  and  1964 


Table  1.  Stations  and  Coordinates  Whose  Results  are  Shown  in  the  PresontReport 
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Station 

Geographic 
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Geomagnetic 

Latitude(°) 
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Ottawa,  Canada 

46N 

75W 

56N 

Slough,  Great  Britain 

51N 

oow 

54N 

Lindau,  Fed  Rep.  Germany 

51N 

10E 

52N 

Moscow,  USSR 

55N 

3?E 

50N 
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For  the  purposes  of  comparing  these  results  with  those  for  other  regions. 
Figure  lb  shows  the  relative  standard  deviation  for  the  same  data  as  shown  in 
Figure  la.  It  can  be  seen  that  the  relative  standard  deviation  is  generally  less  than 
six  percent  implying  that  95  percent,  of  all  observations  lie  within  +  12  percent  of 
the  monthly  median  value. 
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Figure  lb.  Relative  Standard  Deviation  About  the  Observed  Monthly 
Median  foE  at  Ottawa  and  Slough  During  March,  June  and  September 
1958  and  1964 


Figure  2  shows  the  absolute  (2a)  and  relative  (2b)  standard  deviation  about  the 
monthly  median  value  for  foFl  for  June  1958  and  1864  observed  at  Ottawa  and  Slough. 
As  was  the  case  with  foE,  foFl  la  strongly  solar  dependent  permitting  the  calcula¬ 
tion  of  reliable  statistics  only  when  the  solar  zenith  angle  is  relatively  high.  Com¬ 
paring  Figure  2a  with  Figure  la,  it  appears  that  foFl  displays  an  absolute  variabil¬ 
ity  that  to  twice  as  large  as  foE  during  solar  maximum  but  there  is  no  significant 
difference  during  solar  minimum.  The  relative  standard  deviations  (Figures  2b  and 
lb)  indicate  foFl  is  slightly  more  variable  than  foE,  with  the  difference  in  the 
variabilities  of  the  two  regions  being  greatest  during  solar  maximum. 
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Figure  2a.  Absolute  Standard  Devia¬ 
tion  About  the  Observed  Monthly 
Median  foFl  at  Ottawa  and  Slough 
During  June,  1958  and  1964 
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Figure  2b.  Relative  Standard  Devia¬ 
tion  About  the  Observed  Monthly 
Median  foF  .l  at  Ottawa  and  Slough 
During  June.  1958  and  1984 


In  Figure  3  the  local  time  variation  of  the  absolute  standard  deviation  about 
the  monthly  median  F2  region  critical  frequency  is  given  for  three  locations: 
Slough;  Undar.  Germany:  and  Moscow,  USSR.  The  results  are  shown  for  the  four 
months  surrounding  the  solstices  und  equinoxes  for  both  solar  maximum  and  solar 
minimum  conditions.  It  should  bo  noted  that  the  ordinate  in  Figure  3  is  a  factor  of 
10  larger  than  that  in  either  Figure  la  or  2a.  Not  surprisingly,  the  absolute 
standard  deviation  la  larger  at  solar  maximum  (1958)  than  at  solar  minimum 
(1964).  The  results  show  that  at  solar  maximum,  the  variability  of  the  F2  region 
in  generally  larger  for  all  hours  during  equinoxes  than  for  the  time  period 
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Figure  3,  Loca1  Time  Varie'.loa  of  the  Absolute  Standard  Deviation  About  the 
Monthly  Median  foF2  at  Clough,  Lindaw  md  Moaeow  During  March,  June, 
September  and  December,  1058  and  1064 
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surrounding  the  solstices.  During  solar  minimum,  this  doesn't  appear  to  be  the 
case,  however.  In  fact,  the  standard  deviation  appears  to  show  little  seasonal  or 
diurnal  dependence  during  1964,  being  approximately  0.50  MHz  for  all  hours  and 
seasons. 

During  solar  maximum,  the  diurnal  variation  of  the  standard  deviations  at 
Slough,  Lindau  and  Moscow  display  strikingly  similar  behavior  at  all  seasons. 

For  the  equinoctial  months,  the  largest  standard  deviation  (~2.0  MHz)  occurs 
around  0900  h  LMT  and  generally  tends  to  decrease  at  other  hours.  During  the 
northern  solstice,  the  standard  deviation  is  higher  between  the  hours  0000  to  0600 
LMT  than  at  other  times  of  the  day  while  during  the  southern  solstice  the  results 
indicate  the  opposite  behavior. 

To  compare  the  results  of  the  standard  deviations  seen  in  the  F2  region  with 
those  observed  in  the  E  and  FI  regions.  Figure  4  shows  the  relative  standard  de¬ 
viation  for  the  same  times  and  stations  illustrated  in  Figure  6.  It  is  immediately 
obvious  that  the  local  time,  seasonal  and  solar  cycle  dependencies  noted  in  the 
absolute  standard  deviation  art  not  as  apparent  in  the  relative  standard  deviation. 
The  results  in  Figure  4  suggest  that  an  average  value  of  15  percent  provides  a 
good  measure  of  the  standard  deviation  observed  in  the  F2  region  for  all  temporal 
scales.  This  is  two  to  three  times  the  variability  that  is  observed  in  the  E  and 
FI  regions. 

The  fact  that  the  F2  region  is  inherently  more  variable  than  the  E  or  FI 
region  is  certainly  not  a  new  finding  that  is  unique  to  this  study.  However,  the 
evidence  presented  here,  which  is  typical  of  the  mid-latitude  ionosphere,  using 
data  taken  simultaneously  In  time  and  location,  does  provide  some  Indication  as  to 
how  much  more  variable  the  F2  region  is  compared  to  the  ionization  at  lower  alti¬ 
tudes. 

Not  only  are  the  E  and  FI  regions  less  variable  than  the  F2  region,  but  the 
median  value  of  either  foE  and  foFl  at  mid -latitudes  can,  for  the  most  part,  be 
predicted  to  within  an  accuracy  on  the  order  of  five  percent.  (Muggleton,  1972; 
DuCharme  <?t  al,  197 1),  There  appears  to  be  little  doubt,  therefore,  that  propa¬ 
gation  predictions  based  upon  predicted  median  values  of  foE  and  foFl  should  be 
quite  useful  when  applied  to  day-to-day  operational  requirements.  In  the  P’2  region 
the  problem  is  much  more  complex,  The  results  presented  here  show  the  vari¬ 
ability  of  foF2  about  the  observed  monthly  median  value.  Clearly,  median  pre¬ 
dictions  of  foF2  can,  at  least,  lead  to  day-to-day  prediction  errors  on  the  order 
of  0. 8  to  3, 0  MHz  depending  on  local  time  and  phase  of  the  solar  cycle.  Extreme 
caution  should  be  exercised,  therefore,  in  applying  on  n  day-to-day  basis  HF 
propagation  predictions  tliat  depend  on  median  conditions  for  those  radio  trans¬ 
missions  that  are  obviously  controlled  by  the  P2  region. 
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Figure  4,  Local  Time  Variation  of  the  Relative  Standard  Deviation  About  the 
Monthly  Median  r‘oF2  at  Slough,  Lindau  and  Moscow  During  March,  June, 
September  and  December,  1958  and  1984 
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3.  SHORT-TERM  PREDICTION  OF  foF2 

Because  the  F2  region  displays  considerable  variability,  there  have  been  many 
attempts  to  develop  a  short-term  prediction  capability  for  foF2  on  a  time  scale 
shorter  than  a  monthly  median  (Bennett  and  Friedland,  1970).  Most  of  these  studies 
have  been  concerned  with  relating  changes  in  foF2  to  corresponding  changes  in 
selected  geophysical  variables  such  as  the  10. 7  cm  solar  flux  and  the  geomagnetic 
activity  index,  Kp.  The  disadvantage  of  this  type  of  approach  when  applied  to 
timely  predictions  is  that  the  independent  geophysical  variables  upon  which  changes 
in  foF2  depend,  must  themselves  be  predicted.  Thus,  any  errors  in  the  predicted 
values  of  the  geophysical  variables  will  be  reflected  in  the  corresponding  predicted 
values  of  foF2.  A  notable  exception,  however,  is  the  work  of  Gautier  and 
Zacharisen  (1965)  and  Zacharisen  (1965)  who  have  attempted  to  develop  a  short¬ 
term  prediction  scheme  based  upon  statistical  properties  of  the  foF2  distribution. 
Using  auto-correlation  and  cross -corre-  dior.  coefficients  determined  from  devia¬ 
tions  of  foF2  from  the  monthly  mean  value,  they  have  shown  how  much  improvement 
can  be  obtained  using  timely  observation  of  foF2.  In  actual  practice,  however,  the 
monthly  mean  value  of  foF2  is  not  known  in  advance  and  it  must  therefore  be  pre¬ 
dicted. 

In  this  study  we  attempted  to  find  a  short-term  prediction  scheme  for  foF2 
based  upon  immediate  past  observations  of  foF2.  It  was  assumed  that  data  were 
available  at  particular  locations  at  least  24  h  in  advance  of  the  time  a  prediction  is 
to  be  made.  The  problem  then  is  to  find  a  sufficiently  general  prediction  scheme 
that  leads  to  prediction  errors  that  are  smaller  than  those  obtained  using  median 
predictions  over  a  large  temporal  and  geographical  scale. 

To  find  the  optimum  prediction  technique,  various  predictors  were  used  and 
deviations  between  predictor  and  obsarvnt'ou  were  computed  on  an  hourly  basis 
for  each  predictor.  The  standard  deviation  for  each  data  set  was  then  computed 
and  used  to  provide  an  indication  of  the  prediction  error.  The  predictors  used 
were  running  means,  running  medians  and  weighted  means  for  each  hour.  The 
value  computed  from  the  past  data  was  used  to  predict  the  value  of  foF2  for  the 
same  hour  in  3-,  5-,  7-,  9-  and  16 -day  intervals,  A  prediction  was  made  for 
every  hour  of  every  day  and  the  deviations  at  each  hour  were  grouped  into  the 
normal  monthly  intervals. 

Comparing  both  the  standard  deviations  between  predictor  and  observation  as 
well  as  the  number  of  times  a  particular  predictor  came  closest  to  the  observa¬ 
tion  showed  that  the  weighted  mean  predictor  was  the  best  possible  orodi.  'tori  tech¬ 
nique  of  those  that  were  Investigated.  A  simple  weighting  scheme  was  used  in 
computing  tlie  weigliied  mean  predictors  and  was  of  the  form 
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where  m  is  the  number  of  days  used  in  computing  the  weighted  mean  and  D  is 
the  value  of  foF2  the  day  preceding  the  prediction  day;  D_2  is  the  value  of  foF2  two 
days  preceding  the  prediction  day  and  D_m  is  the  value  of  foF2  m  days  preceding 
the  prediction  day.  If  a  value  for  a  particular  day  were  missing,  no  prediction  was 
made  for  that  day. 

Figure  5  shows  the  diurnal  behavior  of  the  standard  deviation  in  MHz  that  re¬ 
sults  using  a  3-,  5  -  and  7  -day  weighted  mean  prediction  for  the  months  indicated 
using  data  observed  at  Slough,  Lindau  and  Moscow  during  1958.  Also  shown  for 
comparison  i3  the  standard  deviation  about  the  observed  monthly  median  value  of 
foF2.  Other  than  for  a  few  hours  centered  near  midday  in  March,  the  standard 
deviation  either  closely  agrees  with  or  is  less  than  the  standard  deviation  about  the 
observed  monthly  median.  Obviously,  if  predicted  monthly  median  values  were 
used,  then  standard  deviation  about  the  predicted  median  would  be  larger  than  those 
computed  about  the  observed  median.  For  all  practical  purposes,  results  shown 
in  Figure  5,  along  with  those  from  other  stations,  indicate  that  using  past  foF2 
data  as  a  predictor,  the  standard  deviation  can  be  reduced  to  or  closely  match  that 
computed  from  the  observed  median.  The  inherent  advantage  in  using  past  observa¬ 
tions  to  predict  future  values  of  foF2  lies  in  the  fact  that  the  parameter  itself  is 
used  in  the  prediction  and  there  are  no  dependencies  on  independent  geophysical 
variables  which  could  be  a  possible  source  of  error. 

Judging  from  Figure  5,  there  appears  to  be  little  difference  in  She  standard 
deviations  using  a  3-,  5-,  or  7  -day  weighted  moan,  and  in  principle,  any  of  the 
throe  should  yield  comparable  results.  In  practice,  however,  it  would  be  desirable 
to  keep  to  a  minimum  the  amount  of  data  needed  to  make  an  accurate  prodiction 
and  it  would  appear  that  the  predictor  requiring  the  smallest  data  bass  would  be 
the  optimum  to  implement.  We  chose  to  uae  the  5 -day  weighted  mean  as  the 
predictor. 

Despite  the  fact  that  observations  of  foF2  can  be  used  to  predict  future  values 
of  ioF2  with  an  accuracy  comparable  to  that  oi'  observed  median,  the  standard  de¬ 
viations,  at  least  during  solar  maximum,  are  quite  large  at  certain  times.  The 
question  arises  as  to  how  useful  are  observations  gathered  on  a  time  scale  less 
than  24  h  in  reducing  the  prediction  error.  To  test  this  out,  it  was  assumed  that 
data  are  available  one  to  six  hours  in  advance  of  the  time  a  prediction  *b  to  be 
made.  These  data  were  used  to  ’update*  the  5 -day  weighted  mean  prediction. 
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Figure  5.  Diurnal  Behavior  at  the  Standard 
Deviation  of  the  Errara  Asaoci&ted  With  a  3' 
5-  and  f-day  Weighted-mean  Prediction  at 
Slough,  f^lndau  and  Moscow  During  1963 
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The  updating  scheme  employed  was  a  straightforward  percentage  change.  If  data 
were  one  hour  old  and  if  the  difference  between  the  observation  and  the  5 -day 
weighted  mean  prediction  was  ten  percent  at  that  hour  then  the  5  -day  weighted  mean 
for  the  next  hour  was  adjusted  by  ten  percent.  Similar  schemes  were  applied  to 
data  that  were  2-,  3-,  .  .  .  6-h  old. 

Figure  6  shows  the  diurnal  variation  of  the  standard  deviation  resulting  from 
updating  the  5-day  weighted  mean  with  data  that  are  1-,  2-.  o-,  4-  and  6-h  old 
using  the  same  stations  illustrated  in  Figure  5,  It  can  be  seen  that  the  amount  of 
improvement  in  the  prediction  (smaller  standard  deviation)  generally  decreases 
with  increasing  time  lag  used  in  the  updating  process.  There  are  times  (0600  h 
LMT  in  March  and  0900  h  LMT  in  December)  when  the  only  improvement  afforded 
is  obtained  using  data  1-h  old  while  at  oth«r  times,  particularly  during  the  midday 
period  at  equinox  and  most  all  hours  during  the  summer  solstice,  data  up  to 
4-h  old,  yield  significant  improvements  in  the  prediction  of  foF2,  During  the 
winter  solstice,  it  appears  that  data  over  2-h  old  applied  in  the  manner  as  des¬ 
cribed  above,  are  of  little  use  in  reducing  the  prediction  error. 

The  effect  of  short-term  data  in  reducing  the  prediction  error  in  foF2  clearly 
has  a  diurnal  dependence.  Data  observed  shortly  before  or  at  local  sunrise  are  of 
little  use  in  updating  the  5 -day  weighted  mean  except  within  1  or  2  h  after  sunrise. 
On  the  other  hand,  foF2  observed  around  0900  h  LMT  can  be  used  to  reduce  the 
prediction  error  for  a  time  period  3  to  6  h  in  advance.  At  periods  surrounding 
local  sunset,  a  behavior  comparable  to  that  at  sunrise  is  seen  .in  the  effectiveness 
of  short-term  observations  of  foF2  to  reduce  the  prediction  error.  Generally,  it 
can  be  said  that  when  the  ionosphere  undergoes  significant  changes  such  as  accom  - 
pony  the  sunrise  and  sunset  transitions,  observations  of  foF2  before  such  transi¬ 
tions  are  of  little  value  when  applied  to  the  updating  scheme  used  here. 

Figure  7  shows  the  comparison  between  the  standard  deviations  obtained  using 
the  media:*.  foF2,  the  5 -day  weighted  mean  foF3  and  the  5 -day  weighted  mean 
updated  using  the  previous  hour's  observation  for  the  1984  data  at  Slough,  Llndau 
and  Moscow.  Generally,  the  8 -day  weighted  mean  yields  a  bettor  approximation 
to  the  observed  daily  hourly  values  than  does  the  observed  monthly  median.  It 
can  be  seen  that  in  many  instances  updating  the  5 -day  weighted  mean  with  data  only 
1-h  old  does  n  A  lead  to  a  significant  reduction  in  the  standard  deviations.  In 
fact,  in  some  instances  (midday'  in  December),  such  updating  causes  an  increase 
in  the  standard  deviation.  The  difference  In  the  effect  of  short-term  updating  in 
foF2  that  is  observed  between  solar  maximum  and  minimum  is  probably  associated 
with  the  fact  that  the  ionosphere  during  soUr  minimum  is  not  subjected  to  a£  many 
of  the  large-scale  perturbations  tlunt  occur  during  solar  maximum  conditions. 
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Flgur«  8.  Diurnal  Variation  of  the  Standard 
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These  perturbations,  due  to  geophysical  phenomena,  such  as  solar  flares  and  geo¬ 
magnetic  storms,  can  lead  to  rapid  changes  in  the  ionosphere  within  a  period  of 
hours.  Thus,  data  obtained  over  five  days  would  be  inherently  more  variable 
during  solar  maximum  than  during  solar  minimum.  This  would  imply  that  data 
taken  in  hourly  increments  should  have  a  larger  effect  on  reducing  the  overall 
ionospheric  variability  during  solar  maximum  conditions  than  during  solar  mini¬ 
mum  conditions. 

It  is  quite  conceivable  also,  that  the  degree  of  improvements  afforded  by  short¬ 
term  updating  is  a  strong  function  of  the  updating  scheme  itself.  We  have  given 
equal  weight  in  the  updating  procedure  to  all  data.  It  is  probably  more  reasonable 
to  weigh  the  effect  of  the  observation  used  in  the  short-term  update  according  to 
the  lag  between  the  time  the  observation  is  made  and  the  time  it  is  applied  to  cor¬ 
rect  the  5 -day  weighted  mean  prediction.  Thus,  smaller  weights  would  be  used  as 
tue  lag  increases.  However,  the  amount  of  improvement  in  the  prediction  by  using 
decreasing  weighting  factors  becomes  correspondingly  small,  thus  diminishing 
the  usefulness  of  short-term  updating. 


4  DISCUSSION 

The  results  described  in  the  previous  sections  clearly  show  that  the  F2  region 
is  far  more  variable  than  either  the  E  or  Ft  regions.  The  relative  variability  of 
ioF2  about  the  observed  median  is  two  to  throe  times  that  of  foK  and  foFl.  Tho 
standard  deviations  of  foE  and  foFl  about  the  observed  monthly  medians  are  such 
♦bet  the  monthly  median  value  of  those  parameters  at  a  particular  hour  provides  a 
reasonable  estimate  of  the  daily  hourly  value  at  mid-latitudes.  This  would  imply 
that  i!F  prediction  iaciuiiquee  such  a3  described  by  Barghausen  et  al,  <1963)  and 
Headrick  et  al,  <19?  1),  which  rely  on  monthly  median  ionospheric  predictions,  can 
be  employed  with  some  confidence  in  day-to-day  operation?  for  those  propagation 
paths  that  are  controlled  by  the  normal  E  and  FI  regions.  For  the  F2  region,  no 
such  statement  can  be  made.  The  standard  deviation  of  foF2  typically  varies  from 
0,6  to  0,7  Mils  during  solar  minimum  conditions  to  0.7  to  2.6  MHz  during  solar 
maximum  Uncertainties  of  this  order  in  foF2  when  applied  to  oblique  transmissions 
via  tho  F2  region  lead  to  uncertainties  In  the  computed  maximum  usable  frequency 
on  the  order  of  1  to  2  Mite  during  solar  minimum  and  2  to  5  MHz  during  solar 
maximum,  assuming  the  monthly  median  value  of  foF2  can  be  accur^ely  predicted. 

To  avoid  errors  associated  with  the  inability  to  predict  the  monthly  median 
value  of  foF2  exactly,  a  prediction  scheme  baaed  on  observations  of  foF2  observed 
over  a  5-day  period  tout  been  shown  to  yield  results  that  ore  comparable  to  or  better 
than  the  observed  monthly  median.  In  addition,  it  has  been  even  that  during  certain 
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times,  particularly  daylight  hours  at  solar  maximum,  observations  of  foF2  ob¬ 
tained  in  hourly,  bi -hourly  or  tri -hourly  intervals  can  be  used  to  produce  signifi¬ 
cant  improvement  in  the  prediction  of  foF2.  Whether  such  a  scheme  is  cost- 
effective,  however,  is  open  to  serious  question.  It  appears  that  using  the  statistical 
type  of  approach  adopted  here,  an  uncertainty  on  the  order  of  0. 5  MHz  will  be  found 
at  all  times.  A  further  limitation  is  that  the  predictions  depend  o.i  observations  at 
particular  locations.  These  predictions  must  then  be  extended  or  extrapolated  to 
regions  where  data  are  not  available,  such  as  over  the  oceans.  This  limitation 
can,  however,  be  minimized  somewhat  by  employing  a  synoptic  mapping  technique 
such  as  described  in  a  previous  report  (Rush,  1972). 

There  appears  to  be  little  doubt,  that  further  work  is  required  in  order  to 
specify,  model  and  predict  the  ionospheric  structure  if  such  efforts  are  to  be  applied 
to  the  prediction  of  day-to-day  HF  propagation  characteristics.  The  most  impor¬ 
tant  area  is  in  the  prediction  of  feF2  itself.  As  means  have  been  devised  that  yield 
predicted  values  of  foF2  with  a  standard  devotion  of  0, 5  to  0,7  MHz  using  simple 
.statistical  techniques,  any  future  approaches,  be  they  mathematical  or  physical, 
must  yield  yet  smaller  standard  deviations  if  they  are  to  be  of  any  further  opera¬ 
tional  use. 
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